Topological floating phase in a spatially anisotropic frustrated Ising model 
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We present new results for the ordering process of a two-dimensional Ising model with anisotropic 
frustrating next-nearest-neighbor interactions. We concentrate on a specific wide temperature and 
parameter region to confirm the existence of two particular phases of the model. The first phase is an 
incommensurate algebraically-ordered (floating) phase emerging at the transition from the param- 
agnetic high-temperature phase. Then the model undergoes a transition to an antiferromagnetically 
ordered second phase with diagonal ferromagnetic stripes (ordering wave vector q = {n/2, 7r/2)). We 
analyze the unconventional features appearing in several observables, e.g., energy, structure factors, 
and correlation functions by means of extensive Monte-Carlo simulations and examine carefully the 
influence of the lattice sizes. For the analytical study of the intermediate phase the Villain-Bak the- 
ory is adapted for the present model. Combining both the numerical and analytical work we present 
the quantitative phase diagram of the model, and, in particular, argue in favor of an intermediate 
topological floating phase. 

PACS numbers: 64.60.De, 64.70.Rh, 75.10.Hk, 05.70.Jk, 75.40.Mg 



I. INTRODUCTION 

One of the most studied classical frustrated spin mod- 
els is the axial-next-nearest-neighbor Ising (ANNNI) 
model. The interplay of the frustrating interactions of 
the model on a square lattice yields new physics in two as- 
pects: new phases emerge and the ordering processes are 
influenced by the competition of the different states. In 
particular Fisher and Selke^ reported on the emergence 
of infinitely many commensurate phases in a certain pa- 
rameter region and Rastelli et al.^ demonstrated similar 
features of the ANNNI model very nicely by applying 
Monte-Carlo (MC) simulations on finite lattices. 

In the present work we examine another classical model 
which has frustrating spin interactions through the di- 
agonal next-nearest-neighbor (nnn) bonds on the square 
lattice. This model was first studied by Fan and Wu.^ 
They found, along the conventional ferromagnetic (FM) 
and antiferromagnetic (AFM) phases, a new phase due 
to frustrating interactions which they called superanti- 
ferroniagnetic (SAF), and it is also more often referred 
as columnar in current literature. This model was ex- 
tensively studied in the recent past^, mainly due to the 
interest to the predicted non-universality of the transi- 
tion into the SAF phase. Even now the model could 
harbor some surprises, as, e.g., the recently found phase 
transition change from the second to the first order at 
strong nnn coupling. We note however that the ear- 
lier work on this model was done almost exclusively for 
the case of equal diagonal couplings. Motivated by some 
real material applications^^ one can generalize the above 
model for the anisotropic case of the diagonal couplings 
of different strength or even of different signs. It turns 
out that the anisotropic nn and nnn Ising model has a 
quite rich phase diagram. In particular, it possesses 
the superferro-antiferromagnetic (SFAF) or 4 x 4 ground 
state phase with ordering wave vector q44 = (7r/2,7r/2). 



After a ^ rotation the ordering pattern of the 4x4 
state becomes equivalent to the antiphase of the ANNNI 
model, as one can see from Fig. 1. From the mean- 
field analysis and Kosterlitz-Thouless-type arguments''''^ 
Chitov and Gros^^ also predicted the incommensurate 
(fioating) phase to be stable within a finite intermedi- 
ate temperature range above the low-temperature 4x4 
phase. The same fioating phase also occurs in the 2D 
ANNNI model. ^"'^ It is characterized by a lack of local 
order parameter and algebraic decay of correlation func- 
tions, modulated by plane wave oscillations with an in- 
commensurate wave vector q depending on couplings and 
temperature. Upon cooling and reaching the boundary 
of the 4x4 phase, the wave vector smoothly evolves 
towards its commensurate value. The incommensurate 
(symmetric) wave vector q^'^ = nn is determined by the 
density of domain walls (n) in the direction of the fer- 
romagnetic diagonals (cf. Fig. 1). The transition from 
the floating to disordered phase occurs via proliferation 
of dislocations (melting) of those walls. It is analogous 
to the Kosterlitz-Thouless transition in the classical XY 
model undergoing through decoupling of topological de- 
fects (vortices). 

The concept of incommensurate phases and phase 
transitions from incommensurate to commensurate or- 
dered states was already discussed in many versions of 
the Ising model in two and three dimensions, see for 
example Refs. 13-16. We should stress two qualitative 
differences between 2D and 3D cases. The 2D float- 
ing phase is characterized by the algebraic order/decay 
of the correlation functions and continuous variation of 
the ordering wave vector with temperature and/or cou- 
plings. There exists a local order parameter in the 3D 
case. Also, as the ANNNI model reveals, the variation 
of the ordering wave vector has discreteness (which sur- 
vives even the thermodynamic limit) referred to as the 
devil's staircase. The number of steps depends on the 
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internal parameters, and thus the staircase reflects the 
intrinsic properties of the system. Such devil's staircases 
have also been observed experimentally (see, e.g., Ref. 17 
and 18). Of course a pure numerical analysis of finite sys- 
tems is often not enough to reveal the true discreteness 
of devil's staircase, and some complimentary analytical 
work is highly desirable. In this context it is worth not- 
ing that the 3D generalization of the present model has 
quite interesting properties. In particular, competing in- 
teractions between stacked Ising planes can results in the 
devil's staircase with respect to the ordering wave vector 
in the stacking direction, observed in the experimental 
work of Ohwada et al.^^ 

Here we will present Monte-Carlo (MC) results of the 
spatially anisotropic Ji-J2-Ising model. In particular 
energies, specific heats, correlation functions and their 
Fourier transform - the structure factor will be discussed. 
Strong numerical evidence for a fioating phase within a 
finite temperature region is given which has not been 
observed before. We show that the Villain-Bak theory^ 
can be straightforwardly extended for the present Ising 
model, giving predictions consistent with the MC results. 

The manuscript is structured as follows: In Sec. II the 
model and its known properties are introduced before we 
present in Sec. Ill the MC and analytical results which 
focus mainly on the description and characterization of 
an intermediate fioating phase. In a concluding Sec. IV 
we discuss our findings in the context of frustrated Ising 
models. 



II. MODEL 

The model is given by summing over all interactions 
between nearest neighbors (nn) and next-nearest neigh- 
bors (nnn) of Ising-spin variables S'i = ±1 on a two- 
dimensional square lattice {N ~ L x L, periodic bound- 
ary conditions): 

H = JiY. S, S, + + ^2 E ^» ^1 ■ (1) 

nn nnn^ nnn^ 

The interaction Ji for nearest neighbors can be chosen 
negative or positive which favors a ferromagnetic or anti- 
ferromagnetic Neel state as configuration of total minimal 
energy, i.e., as a ground state. The J2 are chosen to be 
of opposite sign sgni^J!^) = —sgn{J2) for the two perpen- 
dicular directions. In the following numerical analysis we 
will set always J| = — and will choose IJ2I as energy 
unit, i.e., temperatures and the nn coupling are given 
mostly in units of | J2I. 

Ground States 

For the given parameter set we expect the system to 
order in three different ground states depending on the 
strength of the nn coupling Ji.^^ For large antiferromag- 
netic coupling Ji > | J2 1 the Neel ordered (AFM) state is 




FIG. 1. (Color online) The superferro-antiferromagnetic 
(SFAF) or 4 X 4 ground state of the nn and nnn Ising model 
for -IJ2I < Ji < IJ2I. 



the ground state with a total energy of -Eafm — —2NJi 
with all nnn bonds parallel aligned and all nn bonds an- 
tiparallel aligned. On the other side of the phase di- 
agram a ferromagnetically ordered (FM) ground state 
with E-pM = 2iV Ji for Ji < — IJ2I is realized. In the 
intermediate region the SFAF order yields the lowest en- 
ergy where all nnn bonds are satisfied energetically, i.e., 
antiparallel aligned in one direction and parallel aligned 
in the perpendicular direction, and nn bonds are par- 
allel and anti-parallel aligned alternating in both direc- 
tions, thus, the total energy of the nn sum gives zero 
(see also Fig. 1) and the ground-state energy is given 
by i?sFAF = — 2iV I J2I. The degeneracy of these ground 
states is twofold (FM, AFM) or fourfold (SFAF). At the 
transition points Ji = ±| J2I the ground state is degener- 
ate of order L since every state which is constituted out 
of 4 X 4 plaquettes with a total spin = yields the same 
energy. So, the ground states have long-ranged orders, 
except at the points of quantum criticality. 



Phase Diagram 

The phase diagram for the model was introduced in 
Ref. 12 for varying parameters Ji, Etnd J|. The result 
is a three-dimensional qualitative phase diagram includ- 
ing ferromagnetic and various antiferromagnetic phases 
(Fig. 2 of Ref. 12). In the present work the focus lies on 
a one-dimensional cut through this phase diagram with a 
varying nn coupling Ji and a fixed value = ^J2- The 
finite-temperature phase diagram Tc(Ji/| J2I) was quali- 
tatively sketched in Fig. 6 of Ref. 12. Here we present a 
quantitative diagram in Fig. 2, obtained by the direct MC 
simulations along with analytical work. Note that the 
transition temperatures are invariant under the change 
of the sign Ji —Ji- this was double-checked in par- 
ticular for |Ji| < IJ2I. For the interactions |Ji| > IJ2I 
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FIG. 2. (Color online) Phase diagram of the anisotropic Ising 
model for varying nn interactions J1/IJ2I. Three different 
ground state configurations: Neel order (AFM) for Ji > IJ2I, 
superferro-antiferromagnetic order (SFAF) for — IJ2I < Ji < 
IJ2I and ferromagnetic order (FM) for Ji < — IJ2I are sepa- 
rated by the quantum critical points at Ji = ±| J2I. The three 
phases (disordered, floating, and 4x4) meet at the (exactly 
solvable) Lifshitz point at Ji — 0. Critical temperatures are 
determined using Binder cumulants for |Ji| > IJ2I (energy 
scale is set to | Ji| here) and estimated from energies, specific 
heats and order-parameter behavior for |Ji| < IJ2I (energy 
scale is set to | J2I). The transition temperatures in this regime 
strongly depend on the system size - here L < 200 (see text for 
more details). The red solid line for T^^'^^ is obtained from 
Eq. (10). The dashed blue line for Tf^ is given by the low- 
temperature Eq. (6) and Eq. (11) at 0.782(5) < J1/IJ2I < 1 
and connects Tf^ at J1/IJ2I ^ 0.782(5) with the Lifshitz 
point as a guide to the eye. Note also the agreement with a 
qualitative sketch in Ref. 12. 



the numerical critical temperatures were determined us- 
ing the Binder cumulants for different lattice sizes. ^'^'^^ 
From the intersection point of these cumulants the criti- 
cal point can be estimated. ^^'^^ 

On the other hand for |Ji| < IJ2I, the extraction of 
critical temperatures is more complicated since the order 
parameter for the 4x4 state and its Binder cumulant 
show a strong finite-size dependence. Thus, for the esti- 
mation of the critical temperatures in Fig. 2 energies and 
specific heats were also taken into account. 

The nature of the phase transitions for Ji > | J2 1 and 
Ji < — I J2I is of second order and belongs to the Ising 
universality class. However, according to the prediction 
of Ref. 12, the transition from the high-temperature para- 
magnetic phase to the SFAF state is not direct but rather 
involves an intermediate fioating phase which will be the 
main topic of the following sections. 



III. MONTE-CARLO AND ANALYTICAL 
RESULTS 

The MC simulations allow to calculate the energy, spe- 
cific heat and various order parameters for a wide tem- 
perature range and different parameters Ji for finite lat- 
tice systems. Especially the analysis of the finite-size 
behavior and the extraction of the values for the ther- 
modynamic limit is crucial in the investigation of incom- 
mensurate phases. 

For the simulations we used a Metropolis-single-spin 
update'^'* with an additional exchange MC step.^^"'^^ As 
a starting configuration for several independent runs we 
selected the 4x4 state in the appropriate parameter re- 
gion for all temperatures and performed 10^ thermaliza- 
tion steps. This choice was reasoned in the large energy 
steps between different states in the incommensurate re- 
gion (see below) and already proved itself very successful 
in a similar work on a frustrated Ising model. ^ 



A. Energy and Specific Heat 

The behavior of the energy already indicates that the 
ordering processes for the two phase transitions - to the 
AFM/FM phase and to the 4 x 4 state - differ signif- 
icantly. In Fig. 3 the temperature dependence of the 
energy and specific heat is shown for three cases: For 
Ji = O.2IJ2I, O.8IJ2I transitions to the 4x4 state are 
compared with an Ising-phase transition at Ji = 1.5 | J2I. 

The multiple steps in the energies (Figs. 3(a) and (c)), 
which in addition are shifted for different system sizes L, 
hint towards an unusual ordering in the model that in- 
volves different size-dependent intermediate states. This 
strong size dependence suggests already incommensurate 
ordering. 

Even more prominent than the steps in the energy are 
the peaks in the corresponding specific heat (Figs. 3(b) 
and (d)). They coincide with the steps in the energy 
and prove the existence of several intermediate states. 
The number of these states obviously depends on the 
system size; in particular in Fig. 3(d) for each system 
size (L = 16, 32, 64, 128) more peaks are distinguish- 
able. However, another feature is observed: The posi- 
tion of the last peak, i.e., for the lowest temperature 
(T w 1.23 IJ2I), converges with increasing system size. 
A comparison with the energy curve above (Fig. 3(b)) 
shows that below this temperature the system is ordered 
in the 4x4 state according to the energy value. Thus the 
final transition to the ground state is locked at a finite 
temperature similarly to the conventional transitions into 
the AFM or FM states (compare, e.g.. Fig. 3(f)). 

For a first characterization of the variety of transitions 
we recorded energy histograms. In Fig. 4 we show energy 
distributions at Ji = O.8IJ2I for system sizes L = 64, 
L = 128 and L ~ 256 at different temperatures. As 
expected by the stepwise behavior of the energy, the his- 
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FIG. 3. (Color online) Energies and specific heats for some values of nn coupling Ji > 0. In Figs, (a) and (c) multiple steps 
appear which differ for different system sizes before reaching the ground-state energy -Esfaf ~ — 2|J2|. As a comparison a 
converged energy development for the phase transition to the AFM state is shown in (e) - -Eafm = —2 Ji. The specific heats 
in (b) and (d) show multiple peaks which coincide with the steps in the energies shown above whereas in (f) only a single peak 
is observed which converges at the single transition temperature Tc.afm. 



tograms show double peaks for selected temperatures (at 
the size-dependent transitions temperatures between in- 
termediate states and to the ground state) . The position 
of the ground-state transition is locked for the three his- 
tograms shown in Fig. 4(a) but the peak-to-peak distance 
varies drastically with the systems size: Essentially the 
gap between the prominent peaks in the energy distribu- 
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FIG. 4. (Color online) Energy histograms for some transitions 
at Ji = 0.8 I J2I: Double-peak features are very prominent for 
the ground-state transition (a) and intermediate transitions at 
L = 128 (b) and at L = 256 (c) as well. However, the energy 
gaps are not stable and seem to vanish for all histograms in the 
thermodynamic limit. Note also the different energy scales in 
panels (b) and (c). 



tion shrinks by a factor two while the linear system size is 
doubled. The same behavior can be extracted for the in- 
termediate transitions. In Figs. 4(b) and (c) histograms 
are shown which are recorded at higher temperatures for 
(b) L = 128 and (c) L = 256. The muhiple-peaked 
features are again very prominent but the peak-to-peak 
distances shrink (note the different energy scales for both 
figures) and the temperatures are shifted as well. 

Thus, according to the analysis of energy distribu- 
tions all finite-size transitions show strong first-order 
behavior.^* However, this seems not to hold in the ther- 
modynamic limit since the energy gaps tend towards 
zero. 



B. Order Parameters and Correlation Functions 

To gain further insight into the ordering process of the 
system it is useful to define order parameters and ana- 
lyze their behavior around the phase transitions. The 
order parameters for the FM and AFM phases are read- 
ily defined as magnetization and staggered magnetization 
which can also be expressed via the spin-structure factors 



^(q) 



SiS-j 



(2) 



with ordering wave vectors qpM = (0, 0) and Qafm = 
(tt, tt). The wave vector for the SFAF phase is given by 
q44 = (7r/2, 7r/2); the (square) unit cell has four lattice 
spacings in each direction (see Fig. 1). The square root of 
the normalized structure factor at this wave vector yields 
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FIG. 5. (Color online) Order parameters and Binder cumulants for the 4x4 phase, i.e., moments of the structure factor 
8(^/2, n/2) for the same parameters as in Fig. 3. Although in panels (a)-(d) a strong finite-size dependence is visible the 
convergence of the observables is evident. However, an extraction of the exact transition temperature to the ground state from 
the intersection of the different Binder cumulants is only possible for the AFM case in panel (f). 



a good order parameter. The calculation of this order 
parameter can be implemented also using a staggered 
magnetization 



S{n/2,7r/2) 



N 

/o(r•^r^) 
A /l(r^r«): 



ry)/2]%2 

r^ + l)/2]%2. 



(3) 



(4) 



The modulo operation '%2' yields values zero and one, 
and the two versions /g i account for the degeneracy of 
the ground state, i.e., a shift of all spins by one lattice 
spacing. In addition both states can be flipped com- 
pletely. In Fig. 5 this order parameter and its Binder 
cumulant^°'^-^ 

Ub="^{1-^) (5) 



3(to 



are shown for increasing lattice sizes and nn couplings 
Ji ~ 0.2 I J2I (panels a, b) and Ji ~ 0.8 | J2I (panels c, d). 
As a comparison in panels (e, f) we show the behavior of 
the AFM order parameter and its Binder cumulant; from 
the intersection point the transition temperature can be 
extracted. However, for |Ji| < IJ2I, i.e., for the transi- 
tions to the 4x4 state such an analysis is hampered by 
strong finite-site effects which cause an unusual behav- 
ior at the transition: The cumulants do not intersect in 
a single point and exhibit several dips for intermediate 
temperatures. To properly understand the numerical re- 
sults we must recall that the paramagnetic and the 4x4 
phases are separated by the intermediate floating phase 
which, in particular, does not have a conventional local 
order parameter. 
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FIG. 6. (Color online) Structure factors S{q^,q^) for a series 
of momenta < < 7r/2 vs. temperature. The transitions 
between the intermediate states are clearly visible. Note also 
the stepwise behavior of the signals which is only recognizable 
in the logarithmic scale. 



C. Floating Phase (analytical) 

The key to the analytical treatment of the floating 
phase comes from the observation that the 4x4 phase in 
the ground state has exactly the same pattern (2 rows of 

"up" spins, then 2 rows of "down" spins, or (+ H ) for 

brevity) as the antiphase of the 2D ANNNI model (cf. 
Fig. 1) when viewed along the ferromagnetic diagonals. 
So the Villain-Bak theory,"* developed for that model, can 
be adapted here with minor modifications. In the vicinity 
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of the QCP ( Ji = I J2I) at low temperature our model can 
be mapped onto an effective free-fermionic model which 
accounts for the dynamics of the domain walls. In the 
floating phase these walls are not straight anymore. The 
number density of domain walls k is determined via min- 
imization from the following equation valid at low tem- 
perature: 



s(k)==-2/3(J2- Ji)cxp(2/3J2) 



where 



s{k) 



1 



1 



1 . TTK 

— sm 

TT 1 — K 



(6) 



(7) 



and from the physical meaning of k, it is bound < k < 
1/2. The correlation function in the direction perpendic- 
ular to ferromagnetic diagonals decays algebraically in 
the floating phase 



(5(r' + r)5(r')) - r"" cos(q • r) 



with the power- law index 

■n 



(8) 



(9) 



and the wave vector of oscillations q = ttk determined 
by the density of walls k. This floating phase is bound 
by the low-temperature SFAF phase and the disordered 
(paramagnetic) phase at higher temperatures. The criti- 
cal temperature of the phase transition from the floating 
into the 4x4 phase is given by the Miiller-Hartmann- 
Zittartz method in the framework of the free fermionic 
approximation:"*'^^ 



sinh2^J2 • sinh2^(J2 - Jj) = I , 



(10) 



written here for Ji, J2 > 0. As one can see from Fig. 2, 
this equation agrees nicely with the numerical MC re- 
sults, reproducing even the exactly known Ising value of 
Tc at Ji = (the Lifshitz critical point). Qualitatively, 
the phase transition from the floating into the 4x4 phase 
can be understood as freezing of the domain walls into 

the (-1- H ) structure. 

On the other part of the diagram, with the increase 
of temperature the floating phase becomes unstable and 
undergoes a phase transition into the disordered phase 
when the dislocations start to proliferate in the network 
of the domain walls. In this sense the phase transition 
from the (gapless) floating to the (gapped) disordered 
phase is topological, not accompanied by the appear- 
ance/disappearance of any local order parameter. It is 
analogous to the vortex unbinding transition in the clas- 
sical 2D XY model, and the results within the approach 
due to Villain and Bak can be traced to their coun- 
terparts in the Kosterlitz-Thouless theory. The floating 
phase becomes unstable at the critical value of the walls 
density (or at the wave vector qc — ttKc). It is deter- 
mined from the following ansatz: 



K = 1 - 



1 



V2[l+cxp(-4/3Ji)]i/2 



(11) 
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FIG. 7. (Color online) The critical domain wall density Hie = 
Qc/tt [qc is the critical ordering wave vector) as a function 
of couplings. For each given coupling ratio the wave vector 
grows from q = qc a.t the critical temperature of the transition 
from the disordered to the floating phase T^^ up to g = 7r/2 
at the critical temperature T^^^^ of the transition to the 
4x4 phase. The lower bound of k, fCm = 1 - 1/^ = 0.2929 
and the reentrance boundary (see the text) = 0.3008 are 
also shown. 



The critical temperature Tc and density are found 
from solution of the system of equations (6) and (11). 
The result of the numerical solution of these equations for 
Kf, (and thus for the critical wave vector q^. = ttKc) is given 
in Fig. 7. T^^ of the transition from the floating to the 
disordered phase is shown in Fig. 2 by a dashed blue line. 
The low-temperature equation (6) does not work well at 
J ^ 1^2!, since the curve for Tc "overshoots" the exactly 
solvable Lifshitz point Ji = 0. Contrary to Eq. (10), 
Eq. (6) does not cross over smoothly to the exact result at 
Ji = (note that s(l/2) = —2). This indicates the need 
for a better theory, taking into account, e.g., fermionic 
interactions at arbitrary temperature, which we relegate 
for future work. 

In the vicinity of the QCP the system of equations 
(6) and (11) can be solved analytically. An interesting 
feature of the equations is a small reentrance effect, when 
the floating phase enters slightly into the FM and AFM 
domains |Ji| > IJ2I. The reentrance boundary is given 
by Kc — 1^0- where kq w 0.3008 is the root of s{kq) = 0. 
The critical temperature at the point of reentrance can 
be evaluated as: 



Tc 
J2 



\Jl\ = \J2\ 



-4/ln(2i(Ko-K™)) «1.05 , (12) 



where k,„ = 1 — 1/V2 = 0.2929 is the lower bound of 
K. The bound follows from the stability condition for the 
floating (Kosterlitz-Thouless) phase: 



1] < 1/4 



(13) 



One can also flnd that the critical temperature vanishes 
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at the QCP as: 



where 



^p.2/hif , /, , , 

J2 ' Vkl/J2|-1^' 



A = ^s'(ko)(ko - K„) « 0.04 



(14) 



(15) 



At the QCP K — >■ Km and the correlation function index 
tends to the free fermionic (Ising) value 1] 1/4, while 
in the whole region of the floating phase l/8<^/<l/4. 
The effect of reentrance does not contradict the earlier 
prediction^^ that the floating phase can be adjacent only 
to the SFAF phase (commensurability parameter p = 4) 
of the present nn and nnn Ising model, since there are 
no commensurate- incommensurate transitions with small 
< 8, because then 77 = would violate the condi- 
tion (13).'^'''''® Numerically, the reentrance is very small, 
and the reentrant regions do not overlap with the FM 
or AMF phases. It could well be an artefact of the low- 
temperature approximation (6), however, for a reliable 
test by MC simulations parameter region is too small. 



D. Floating Phase (Monte-Carlo) 
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FIG. 8. (Color online) Results for the fit parameters of Eqn. 8 
when applied to the finite systems for L — 32 at Ji — —0.5 | J2I 
and for L = 64 at Ji = -f0.8|J2|. In the upper panel the 
development of the wave number is shown which reflects the 
behavior of the structure factors given in the left-hand side 
of Fig. 9. In the lower panel the algebraic decay exponent is 
shown: For the intermediate phase a saturation at different 
levels can be observed before it holds 77 = in the ground 
state. (Values in the vicinity of transitions are left out since 
there the fitting fails.) 



To verify the predictions for the floating phase, we cal- 
culated the structure factors for symmetric momenta on 
the line qfm = (0,0) — )• q44 — > Qafm = We ob- 

served finite signals for some structure factors depending 
on the sign and value of Ji and more important also de- 
pending on the lattice size. Exemplary we show several 
structure factors at Ji = — O.8IJ2I in Fig- 6 in a loga- 
rithmic scale for a broad temperature range. Due to the 
computational effort here only L = 64 is chosen, but nev- 
ertheless a cascade of signals can be seen which starts at a 
critical momentum < 7r/2 and ends up in the ground- 
state order parameter with qf^ — tt/2. As we show above, 
the Villain-Bak theory for this model predicts that the 
intermediate floating phase is characterized by a single 
wave vector q while in the MC results for all momenta 



(a"" ) 



and 



Qc < qfnc. < 

7r/2 < gf„,. < q^^ 



for Ji < 
for Ji > 



(16) 



a non-zero contribution can be observed. The left-hand 
side of Fig. 9 shows the structure factors for several q's 
also in a color-coded plot for L = 32 at Ji ~ — 0.5| J2I 
(upper left panel) and i = 64 at Ji = +0.8| J2I (lower left 
panel). The behavior described in Eqn. (16) can be seen 
very nicely in these plots. In conclusion in the interme- 
diate phase the structure factor shows finite signals for 
certain momenta which depend (i) on the sign and magni- 
tude of the nn interaction and (ii) on the lattice size. This 
can be explained by the simple fact that for a discrete 
lattice the number of moments qi^c. between qdJi) and 



944 = 7r/2 increases with the system size. Furthermore 
in Fig. 6, the stepwise behavior of the structure factor at 
the transitions between different momenta is very promi- 
nent. However, this feature is also due to the discrete 
spectrum of the momenta on a finite lattice. As already 
observed in the energy histogram finite-size effects play a 
crucial role in the intermediate phase. In the thermody- 
namic limit the temperature dependent momentum q(T) 
locks smoothly into the ground-state value q44. 

The size dependence of the signals in the structure fac- 
tors show already strong indications for unconventional 
ordering in the system before the actual ground state is 
reached. A further analysis relies on the evaluation of 
correlation functions. On the right-hand side of Fig. 9 
we present correlation functions calculated in one row or 
column of the lattice for L 32 at Ji = — 0.5| J2I (upper 
right panel) and L = 64 at Ji = +O.8IJ2I (lower right 
panel). These correlation functions clearly show non- 
fitting oscillatory behavior in the region of higher tem- 
peratures before the system orders in the ground state 
with its four-site period. An analysis of these oscilla- 
tions by fitting Eqn. (8) at each temperature step yields 
a good agreement in a wide temperature range. The re- 
sulting decay exponent rj and the wave number q can be 
extracted. Since the correlations are calculated in one 
dimension only, the wave vector is reduced to one com- 
ponent here. Furthermore, this wave number obviously 
coincides with the symmetric entries of the wave vector 
of the corresponding structure factor which are shown 
on the left-hand side of Fig. 9. The results for q (upper 
panel) and rj (lower panel) are given in Fig. 8. 
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m %/2 % 3/471 S(q) 5 10 15 20 25 303.3. 

q" (= q^) distance r 

(b) L = 64; Ji = +O.8IJ2I 



FIG. 9. (Color online) Structure factors (left) and correlation functions (right) calculated from two different systems shown 
in color code. The stepwise behavior of the structure factors coincides with the distinguished temperature regions in the 
correlation functions, and the features clearly depend on Ji and L. The critical momentum qc is determined numerically as 
described above and shown in Fig. 7. 



In the plot fitting parameters in the vicinity of tran- 
sitions and for high temperatures are left out since at 
these points the fit fails and does not yield meaningful 
results. The upper panel (wave numbers) only repro- 
duces the results from the structure factors but the lower 
panels shows that in the intermediate phase the behav- 
ior of a floating phase is recovered. The fits of the decay 
exponent give reasonable (taking into account the strong 
finite size effects) results in the floating phase. At low 
temperatures the exponent saturates at 77 '--^ 0.2, while 
below Tc into the 4x4 phase with the long-range order 
parameter it adopts the value 77 = 0. 



IV. CONCLUSION 

An anisotropic version of the frustrated J1-J2 Ising 
model was investigated using mainly MC simulations, 
supplemented by analytical treatments. In particular the 



phase transition into an antiferromagnetic ground state 
constituted of two sublattices in collinear order (4x4 
phase) was analyzed. It was predicted earlier by Chi- 
tov and Gros^^ along with an unconventional floating 
phase which appears for intermediate temperatures be- 
fore the ground state is reached. For analysis of the 
floating phase the Villain-Bak theory^^^ with some small 
modifications was utilized for the present model. In the 
floating phase the theory predicts no local order parame- 
ter, algebraically decaying correlation functions with the 
power-law exponent rj, modulated by a plane wave with 
a (single) incommensurate wave vector q depending on 
couplings and temperature. The theory also allows to 
quantitatively describe the smooth evolution of the wave 
vector from its critical value Qc at the critical temperature 
between the disordered and floating phases towards its 
commensurate value at the boundary of the 4x4 phase. 
The critical temperatures between these phases are also 
evaluated. Qualitatively, the theory gives the picture of 
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the transition from the floating to disordered phase via 
proliferation of dislocations of the domain walls, which 
is analogous to the Kosterlitz-Thouless transition of vor- 
tices in the classical XY model, while the transition into 
the commensurate phase occurs via freezing of the do- 
main walls into the 4x4 structure. 

The nature of the intermediate phase was also ana- 
lyzed using correlation functions and the corresponding 
structure factors in the MC simulations. It appears in 
the simulations of the finite-size system that the phase is 
not described by a single wave vector but rather a set of 
neighboring wave vectors. The transitions in finite lat- 
tices between states with different momenta are sharp, 
however, in the thermodynamic limit energy gaps seem 
to vanish. Thus the phase is best described by a temper- 
ature dependent wave vector whose discrete spectrum is 
smeared out in the thermodynamic limit. The momen- 
tum varies from a starting vector (\c{Ji/\J2\) and locks 
smoothly into (7r/2,7r/2). The nature of this phase was 
further analyzed by fitting the correlation functions by a 
combination of algebraic decay and oscillatory behavior. 
The agreement in the intermediate phase is very good 
and we conclude from this result that the state is best 
described by a floating phase consistently with the ana- 
lytical predictions. The model's phase diagram in Fig. 2 
summarizes most of our MC and analytical results. 

As for the further work, we expect a very interest- 
ing behavior of the present model when the transverse 



fleld is included, especially near quantum criticality.^^ 
Another interesting direction is the 3D generalization of 
the model, where the devil's staircase, similar to the one 
in the 3D ANNNI model, ^'^ is expected. It appears that 
the devil's staricase in the 3D extension of the present 
model was already observed experimentally.^^ 
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